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Комплексний аналіз динаміки 
сонячної активності 


Проведений комплексний аналіз часового ряду усереднених місячних значень сонячної активності, 
який включає гармонічний аналіз, К/5-аналіз, аналіз методами хаотичної динаміки та оптимізоване 
фільтрування з використанням дискретного перетворення Фур'є. В результаті була виділена детермінована 
складова ряду у вигляді полігармонічної моделі. 


Вступ 


У зв'язку із глобальним потеплінням дослідження змін клімату та спроби їх 
прогнозування є найважливішою задачею, яка стоїть перед людством. Одним з головних 
факторів, які впливають на систематичні зміни погоди і клімату, є зміни сонячної актив- 
ності (1). Наша робота присвячена аналізу динаміки сонячної активності, вираженої 
через ряд чисел Вольфа за період 1750 - 2006 рр. |2). Цей ряд містить 3084 числа 
(257 років х 12 місяців), які є усередненими місячними значеннями сонячної активності. 
Найважливішим і поки що невирішеним питанням є питання про характер спостере- 
жуваної динаміки сонячної активності. Перший варіант - детермінована поведінка в 
рамках нелінійної (хаотичної) динаміки. Другий варіант - стохастична поведінка згідно з 
моделями теорії випадкових процесів. Недавні дослідження, проведені в рамках неліній- 
них динамічних систем (3, |41, дають ряд доказів присутності низькорозмірної хаотичної 
поведінки системи сонячної активності. Проте не можна вважати ці докази остаточними, 
оскільки немає універсального методу, який би надійно розділяв частотний спектр на 
кольоровий шум і приховані динамічні процеси. Нами запропонована нова методика 
дослідження часового ряду сонячної активності, яка включає фільтрування з викорис- 
танням дискретного перетворення Фур'є у поєднанні з В/5-аналізом та гармонічним 
аналізом часового ряду. 


К/5-аналіз часових рядів 


Одним із потужних методів дослідження довготермінової пам'яті у часових рядах 
та з'ясування ступеня їх детермінованості є В/5-аналіз (КК - розмах, 5 - стандартне 
відхилення). Метод не використовує жодних гіпотез про вигляд функції розподілу. 
Запропонований Херстом алгоритм В/З-аналізу |5| обчислює середнє значення 
показника Херста Н є (0,1). Якщо показник Н приблизно дорівнює 0,5 - це свідчить про 


те, що ряд описує випадкове блукання. Якщо 0,53 Н «1 - ряд є додатнокорельованим, 
або ж трендостійким. Чим більшим є показник Херста, тим кращою є трендостійкість 
часового ряду і тим більше підстав для його надійного прогнозування. Якщо виходити з 
гіпотези про адитивну модель ряду у вигляді «детермінована складова» "| «шум», то 
буде логічним використовувати показник Херста як індикатор співвідношення 
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«сигнал-шум». Проведені нами розрахунки показника Херста для ряду середньомісячних 
даних сонячної активності дали значення Н - 0,77. Таке значення Н є доказом присут- 


ності ефекту довготермінової пам'яті і помітного ступеня детермінованості системи. 


Гармонічний аналіз ряду сонячної активності 


Циклічність може служити ознакою прихованої динаміки процесу. У випадку 
сонячної активності загальновизнаним є 11-річний цикл. Періодограма, побудована в 
середовищі 5ТАТІЗТІСА, вказує також на наявність циклів з періодами 10 1 100 років. 
Комбінований метод найменших квадратів та повного перебору |6| дозволяє побудувати 
набагато більш точну періодограму ряду середньомісячних значень сонячної активності. 
Цей метод дозволяє послідовно виділити цикли, які є найбільш характерними для 
даної динамічної системи (гармоніки з найбільшими амплітудами). У загальному вигляді 
феноменологічна полігармонічна модель часового ряду сонячної активності має 
вигляд (па - кількість гармонік) 


т 
х зі з Уа, зш( -к Ф,). (1) 
ізі 
Тут РИ - модельне значення сонячної активності, Ф, - середнє значення сонячної 
активності за період спостережень, а, - амплітуда 1-ої гармоніки, (2, - циклічна 
частота, Ф; - фаза коливань, / - час. Значення параметрів 1-ої гармоніки (1 - |, 2, 3,...) 
послідовно визначалися з умови мінімуму функціоналу похибки моделі 


Т тах 
- 2 - 
у з У (х, за; "віщо н Ф,) - Б.) -» тій. (2) 
те0 
Тут Х, - фактичні значення сонячної активності; х, - Х, - х,, х'", є сумою попередніх 
і - 1 гармонік. Зауважимо, що ф, - 0 при 121, але, внаслідок похибки обчислювальної 


схеми, ці параметри можуть приймати значення, які дещо відрізняються від нуля. 
Результати гармонічного аналізу динаміки сонячної активності наведені у табл. І. 
Аналіз амплітуд гармонік вказує на те, що основними є перші три гармоніки. 


Таблиця 1 - Гармонічний аналіз ряду середньомісячних значень сонячної активності 


Період (роки) |Початкова фаза (роки)| Амплітуда | Вільний член 
1-а гармоніка 11,0 1,2 29,11 52,43 
2-а гармоніка 10,0 0,4 24,59 0,09 
3-я гармоніка 10,6 9,0 18,39 0,05 
4-а гармоніка 103,8 10,8 15,54 -1,33 
5-а гармоніка 11,8 5.3 15,46 -0,15 


Виникає питання - наскільки адекватно відображають виділені гармоніки 
приховану динаміку системи? Якщо виділені гармоніки є проявом внутрісистемної 
динаміки, вилучення їх повинно робити ряд більш хаотичним, тобто зменшувати 
значення показника Херста. В результаті виконаного В/5-аналізу нами отримані 
наступні значення показника Херста: початковий ряд: Н «0,77; ряд перших 
залишків, утворений після вилучення першої гармоніки: Н «0,82; ряд других 
залишків: Н «0,85. Оскільки виділення перших гармонік не зменшує показник 
Херста, то немає підстав твердити, що полігармонічна модель (1) з параметрами табл. 1 
адекватно моделює детерміновану складову ряду. Причина може полягати у значній 
зашумленості початкового ряду. 
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Згладжування часових рядів 


Кожен часовий ряд, який генерується природною чи економічною системою, є 
значною мірою «зашумленим» стохастичними збуреннями. Для виявлення прихованої 
динаміки слід перш за все згладити часовий ряд, тобто відфільтрувати високочастотні 
шуми. Ефективне фільтрування підсилює основний «сигнал» і приглушує «шум». Для 
досягнення такого ефекту необхідно, щоб фільтр був узгодженим із сигналом |7|. Це 
означає, що його ефективна ширина повинна бути близькою до ширини піків та впадин 
сигналу. Оскільки ми не знаємо форми справжнього сигналу, то будемо підбирати 
ширину згладжуючого вікна з деяких додаткових міркувань. По-перше, значення 
ширини повинно перевищувати 12 (кількість місяців у році). По-друге, відфільтро- 
ваний ряд повинен володіти властивістю зменшення коефіцієнта Херста після виділення 
перших гармонік. По-третє, якщо ми хочемо отримати докази присутності динамічного 
хаосу, необхідно виконання умови насичення кореляційної розмірності. 

Одним з найбільш поширених методів фільтрування сигналу є метод віконного 
перетворення Фур'є |3). Для покращення частотної характеристики фільтра використо- 
вують зважуючі коефіцієнти, які плавно зменшують до нуля досліджуваний сигнал на 
краях вікна. Якщо вважати, що сторонні шуми є гауссовим випадковим процесом, то опти- 
мальною формою зважуючого вікна буде гауссіан. Ядром згортки виступає функція 

2 


-х 
К(ху) з з (3) 


Тут х - елементи ряду сонячної активності, с - параметр, що визначає ширину 
вікна. Чим більша ширина гауссової кривої, тим сильнішим буде ефект згладжуван- 
ня. Повна ширина зважуючого вікна визначається співвідношенням |З8| 


у, «ду е2с--/815(2). (4) 


Алгоритм фільтрування складається з наступних кроків: 
. Перетворення Фур'є початкового ряду. 

. Перетворення Фур'є ядра згортки. 

. Множення отриманих перетворень. 

. Обернене Фур'є-перетворення добутку. 
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Рисунок 1 - Середньомісячні кількості сонячних плям 
(початкові дані - точки; згладжені дані ( у/ - 33) - товста лінія) 
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Розрахунки були виконані нами з використанням пакету Мафетайса 4.1. Ми 
вибрали параметр ширини вікна у - 33, що вдвічі перевищує значення, вибране 
авторами роботи |3)|. Це значення було визначене нами як оптимальне після ряду 
комп'ютерних експериментів. Наш вибір буде обгрунтований в наступних вик- 
ладках. Зараз відмітимо лише те, що повна ширина зважуючого вікна (66 точок) 
дорівнює половині довжини основного сонячного циклу (11912 - 132). Вигляд 
згладженого ряду представлений на рис. І. 


Побудова фазового портрета системи 


Наступним етапом нашої роботи є дослідження відфільтрованого ряду з метою 
виявлення поведінки, характерної для хаотичної динаміки. Для цього необхідно виконати 
побудову фазового портрета системи та обчислити кореляційну розмірність а. Першим 
кроком у дослідженні динамічних характеристик скалярного часового ряду є реконструк- 
ція І) -вимірного фазового простору «методом затримок» |91. Необхідно встановити 
мінімальну розмірність простору вкладення /), яка забезпечить топологічно адекватне 
відображення фазового портрета системи (фазова траєкторія не повинна мати самопере- 
тинів). Для визначення розмірності вкладення використовується метод «фальшивих 
сусідів» (10). Необхідно збільшувати розмірність фазового простору до тих пір, поки 
середня відстань між векторами - найближчими сусідами - перестане змінюватися. 
Результати наших досліджень розмірності вкладення представлені на рис. 2. Значення 
параметра затримки г «- 30 є близьким до значень, використаних в роботі (31. 
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Рисунок 2 - Залежність відносної зміни середньої відстані між сусідніми точками К 
від розмірності фазового простору Р. Час затримки (в місяцях) г - 30 


Як видно з рисунка, насичення наступає при 0 « 4 - 5. Проте, враховуючи той 
факт, що згладжування не усуває всі шумові ефекти, можна стверджувати, що 
фактичне значення розмірності вкладення 2) - 3. Аналогічні результати отримані в 
роботі |3| при значенні згладжуючого параметра у» «15. 
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Рисунок 3 - Проекція фазового портрета, отриманого з використанням параметрів 
у «33, те30 


Фазовий портрет, який відповідає параметрам у - 33, т - 30, зображений на 
рис. 3. Згладжування з вилученням високочастотних флуктуацій допомагає відновити 
приховану динаміку системи. Не виключено, звичайно, що згладжування може вида- 
лити частину оригінальної динаміки, але зрозуміло, що ця процедура не може привнести 
елементи нелінійної динаміки у згладжені дані. 


Комплексний аналіз згладженого часового ряду 


Виконаємо комплексний аналіз згладженого ряду середньомісячних значень 
сонячної активності, яке включає гармонічний аналіз та В/З-аналіз. Результати 
гармонічного аналізу наведені у табл. 2. 


Таблиця 2 - Гармонічний аналіз згладженого ряду ( у - 33 ) середньомісячних 
значень сонячної активності 


Період (роки) | Початкова фаза (роки)| Амплітуда |Вільний член 
1-а гармоніка 11,0 1,2 17,36 52,43 
2-а гармоніка 104,0 10,8 15,15 0,09 
3-я гармоніка 10,6 9,0 13,62 0,05 
4-а гармоніка 53,1 36,7 11,42 -1,33 
5-а гармоніка 10,1 9,2 10,41 -0,15 


Результати виконаного нами КВ/5-аналізу є наступні. Початковий ряд: 
Н «0,82; ряд перших залишків: Н «0,76; ряд других залишків: Н «0,66; ряд 
третіх залишків: Н - 0,66. Як бачимо, згладжування підвищило значення показника 
Херста (від 0,77 до 0,82), що свідчить про ефективне видалення високочастотних 
флуктуацій. Виділення перших двох гармонік із згладженого ряду монотонно знижує 
показник Херста до 0,66, що є свідченням адекватності полігармонічної моделі для 
згладжених даних. Все це підтверджує правильний вибір параметрів нашої моделі. 
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Отже, для ефективного моделювання та прогнозування динаміки сонячної 
активності необхідно застосувати полігармонічну модель виду (1) із параметрами, 
наведеними у табл. 2. Крім 11-річного циклу найбільш помітними циклами сонячної 
активності є цикли з періодами 104 та 10,6 року. Щодо останнього циклу, то він може 
бути результатом накладання двох гармонік з періодами 11,0 років та 10,1 року. 

Застосування комбінації гармонічного та В/5-аналізу до згладженого ряду 
середньомісячних значень сонячної активності, використаного в роботі |3| (у  - 15), 
дало наступні результати. Початковий ряд: Н - 0,73; ряд перших залишків: Н « 0,79; 
ряд других залишків: Н «0,81; ряд третіх залишків: Н - 0,93. Зростання показника 
Херста свідчить про недостатнє фільтрування шумів при півширині вікна у - 15. 


Визначення кореляційної розмірності 


Для випадкових процесів залежність кореляційної розмірності від розмірності 
простору вкладення 0) є монотонно зростаючою (11). Насичення залежності «4(Д) є 


першою ознакою динамічної поведінки системи. Для визначення а зазвичай використо- 
вують процедуру Грасебергера -- Прокаччі (12). В основі процедури визначення кореляцій- 
ної розмірності лежить поняття оче ді інтеграла С, (т): 


СЯ і Ук РВЕ (5) 
і,/-І 
Тут О(7) - функція Хевісайда: О(г) - 1, якщо г» 0; 0(/) - 0, якщо " «0. При малих 
г кореляційний інтеграл має вигляд |10)) 
Св Аг? (6) 
або, після логарифмування, 
ішСь(гу- Аза. (7) 
Згідно із рівнянням (7) на графіку залежності Іп(С),) від їй г повинна бути прямолінійна 
ділянка, яка є свідченням самоподібної геометрії фазової траєкторії системи. Будемо 
збільшувати розмірність простору вкладення Д). Починаючи з деякого значення 0) 
показник степеня а! перестає змінюватися (в межах похибки обчислень). Тоді остаточне 
значення кореляційної розмірності а визначають із графіка як середнє значення 
нахилу прямолінійної ділянки графіка, який відповідає даній розмірності вкладення 
Р, при якій наступає насичення. Похибка визначається як різниця між максимальним та 
мінімальним значеннями нахилу прямолінійної ділянки. 

Застосувавши процедуру Грассбергера - Прокаччі до згладженого ряду се- 
редньомісячних значень сонячної активності (у -33), ми отримали значення 
кореляційної розмірності а - 2,16. 0,05. Насичення залежності 4()) наступає вже 
при розмірності вкладення 10) «3-4 і це є свідченням динамічного характеру 
системи. Зауважимо, що для добре вивченої системи Лоренца П 3, розмірність 
вкладення якої дорівнює трьом, значення кореляційної розмірності становить 
д-«2,06. Тому логічним буде припустити, що розмірність вкладення для системи 
сонячної активності також дорівнює трьом, що підтверджує висновки, отримані 
вище методом «фальшивих сусідів», 1 співпадає з результатами роботи |З. 


Висновки 


Нами виконаний комплексний аналіз динаміки сонячної активності, відображеної 
часовим рядом чисел Вольфа. Розроблена оптимізована процедура згладжування 
часового ряду, яка ефективно видаляє високочастотні шуми. Оптимальна ширина 
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згладжуючого вікна визначається за результатами послідовного застосування гармонічного 
аналізу та В/5-аналізу до початкового часового ряду. Запропонована феноменологічна 
полігармонічна модель сонячної активності, адекватність якої підтверджена результатами 
К/5-аналізу. 

Дослідження згладженого ряду методами хаотичної динаміки привели до 
наступних висновків: 

1. Швидке насичення залежності «())) свідчить про динамічний характер 


системи. 
2. Розмірність вкладення системи 0 «3. 
3. Кореляційна розмірність а «2,416 60,05. 
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П.М. Грицюк 

Комплексньй анализ динамики солнечной активности 

Проведен комплексньгюй анализ временного ряда усредненньх месячньх значений солнечной активности, 
которьій включаєт гармонический анализ, В/5-анализ, анализ методами хаотической динамики и 
оптимизированное фильтрование ряда с использованиєм дискретного преобразования Фурье. В результате 
бьтла вьіделена детерминированная составляющая ряда в виде полигармонической модели. 


Реїо М. Нтуізунк 

Сотріех Апаїузіз ої 5о0Іаг Асбїуїку Бупатіся 

Тре соппріех апаЇузія ої їбе те 5егіе5 ої Пе 50Іаг асіїутку плопібіу ауегаре уашез іпстпдтпо, їпе Багтопіса! 
апаіузіє, В/5 апаїузіз, їпе апаЇузі5 бу сраобіс дупатіся пекфодз апа бе орбїтігей Негіпе млі ц5еад ді5сгеїс 
Еоигіег (гап5їогт ухаз саттіва оці. Аз а гезиії Бе те 5егіе5 деїегтаїтед сотропегі а5 а роїубагтопіса! плодєї 
ууа8 ехітастіва. 


Стаття надійшла до редакції 23.11.2007. 
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